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Abstract 



For = 1, 2, . . ., let be a simple random sample of size n = 
from a population An of size A^, where < n < A^. Then with 
/at = n/N, the sampling fraction, and 1a the inclusion indicator that 
A G Sn, for any H C An of size A; > 0, the high order correlations 



depend only on k, and if the sampling fraction /tv — )• / as — t- oo, 
then 



where Z is a standard normal random variable. This proves a conjec- 
ture given in [2]. 
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1 Introduction 



Simple random sampling is doubtless one of the most often used tools in 
statistics [6], and it might appear that nothing new regarding this simple 
scheme remains unexplored. With < n < A^, by a simple random sample 
of size n from a set of size we mean the random subset oi Sn oi An 
with distribution 

P{SM = r) = \ (nf rcAN,\r\=n 
{ U otherwise, 

where |r| denotes the size of the set r. It is easy to see that all individuals 
in An have an equal chance of being included in the sample, and that in 
particular, for A G An, 

Tl 

El A = In where In = j^, (2) 

where the inclusion indicator 1^ takes the value 1 when A & Sn and the value 
otherwise. The value Jn is known as the sampling fraction. Likewise, from 
(II]) one can show directly that the inclusion indicators 1a and 1b, for distinct 
individuals A and B in An, are negatively correlated, as the inclusion of A 
leaves less room in the remaining sample for B. That is, 

E1a1b= ''^;-^\<(-Y = E1aE1b, 
^ ^ N{N-l) \nJ ^ ^' 

or, considering the two way correlation 

Corr(lA, 1b) = E {{1a - Jn^Ib - In) , (3) 

we have 

/-i t \ -n{N-n) 
Corr(U, 1b) = ^,^^_^y (4) 

However, the higher order correlations of simple random sampling, which 
arise in some applications [2] and exhibit rather interesting behavior, are 
virtually unknown. 

To consider such correlations of higher order, generalizing ([3]), for any 
H C An of size \H\ = k,0 < k < n, we define 

Corr{H) = E(ll{lA-fN)]- (5) 
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We see from ([T]) that the probabihty that all individuals in H are included 
in the sample is 

XAeH / \n) 

which we note only depends on fc, and not on which individuals comprise the 
set H . Hence, arguing either directly using (jH]), or by noting the indicators 
{lyi, A G An} are exchangeable, we conclude that Corr(if) depends only on 
the size k of if, and hence we denote it by Corr(A;). 

In [2] the high order correlation of rejective sampling was studied in order 
to determine the asymptotic properties of a generalized logistic estimator. 
In rejective sampling [S], each individual A in An is associated with a non- 
negative weight wa, and the probability of sampling a set r C of size n 
is given by 

P{Sn = r) = — — where Ws = ]^Wj. 

We note that the high order correlations of rejective sampling may be defined 
exactly as in ([5]), and that simple random sampling is the special case of 
rejective sampling, taking all weights equal. 

Critical in the asymptotic analysis in [2] was the fact that under some 
stability conditions on the weights, the second and third order correlations 
of rejective sampling decay at rates A^~^ and A^~^, respectively, that is, that 

lim A^Corr(2) = 0(1), and lim N^Corr{3) = 0(1). 

Checking for the special case of simple random sampling when /tv — )■ / as 
— 7- oo, equality @ implies A^Corr(2) — )■ /(I — /), and further direct cal- 
culation for correlations up to order 9 obtained by expanding the expression 
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in definition and using yields 



A^Corr(2) 




/(/-I) 


A^^Corr(3) 


— > 


2/(/-l)(2/-l) 


A^2Corr(4) 




3f(/-lf 


A^3Corr(5) 




20f (/- 1)^(2/ -1) 


A^3Corr(6) 




15f (/-If 


A^^Corr(7) 




210/3(/- 1)3(2/ -1) 


A^^Corr(8) 




105/^(/ - 1)^ 


A^^Corr(9) 




2520/^(/ - 1)^(2/ - 1) 



Perhaps the most surprising feature of these correlations is that their rate of 
decay depends on the parity of the correlation order, in particular, one can 
easily conjecture that 



Theorem 4.1 of [2] shows that ([7]) holds quite generally for rejective sam- 
pling, and therefore for simple random sampling in particular. Application 
of this theorem sufficed to complete the asymptotic analysis required in [2] 
for rejective sampling. 

Another feature of the simple random sampling correlations is also easy 
to conjecture, that their scaled limits are equal to a constant depending on 
k, times the factor f{f — 1) raised to {k — kmod2)/2, and if k is odd, times 
the additional factor (2/ — 1). Hence, one need only determine the constants 
to completely specify their limiting behavior. On the basis of the above 
observations and the constants corresponding to the even and odd values of 
k up to 9, that is, the sequences 1, 3, 15, 105 and 2, 20, 210, 2520, respectively, 
a conjecture was put forth in [2], which is now validated by the following 
theorem which is proven in this analysis. 

Theorem 1.1. For N = 1,2,... let Sn be a sequence of simple random 
samples from populations An of size N, whose sampling fractions f^ obey 



N{k+kmod2)/2^^^^^^^ = 0(1) A; = 1, 2, ... . 



(7) 



fN^fe (0, 1) 



as iV — )■ oo. 



Then 



lim 




l^']k/2-^2k^ A; even 
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and 

lim Ar('=+i)/2Corr(fc) = [/(/ - l)]^^-^'^/^i2f - - l)EZ^+\ k odd 

N^oo 3 

where Z is a standard normal variate. 

We recall that the standard normal variable Z is the one with distribution 
function ^ 

P{Z <x) = [ e-"'/2^M, 

v27r J-oo 

or, equivalently, moment generating function M{t) = Ee*^ given by 

M{t) = 

Substituting jl for u in the expansion e" = Xljlo'^Vi'' fi'^d 

and hence, as EZ^ = M'^'^^(O), we conclude that Z has vanishing moments 
of odd order, while 

EZ^ = , , . = TT i = -rr, — — ^ for k even. 

2^/2 W2 ! k(k-2)---2 

^ ' ' 0<j<k,jodd ^ ' 

The appearance of the moments of the standard normal variable in count- 
ing contexts is well known. In particular, in the On-Line Encyclopedia of In- 
teger Sequences [7] there are roughly twenty combinatorial structures listed 
which are counted by the sequence EZ'',k = 0,2,4,..., such as the num- 
ber of perfect matchings in the complete graph K{2n), and the number of 
fixed-point-free involutions in the symmetric group S'2„. However, in con- 
trast, the sequence (l/3)(/c — 1)EZ^~^^, k = 1, 3, 5, . . . has only two listings, 
and of a more abstract nature, in particular, the Ramanujan polynomials 
~V^n+2(^ + 2,x) evaluated at 1, and, with offset 2, the second Eulerian trans- 
form of 0, 1, 2, 3, 4, . . .. In addition to providing a concrete situation in which 
this lesser known second sequence appears, the present work also demon- 
strates that there exists a connection between it and the one much better 
known. At the same time, we remark that our proof of Theorem 1 1.1 1 is purely 
computational, and that it would be enlightening to also complete the argu- 
ment in a manner which makes the appearance of the normal moments more 
natural. 
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2 Main Result 



The proof of Theorem 11.11 requires a few identities for the Stirhng and 
Bernoulh numbers, which can be found in Letting {x)j denote the faUing 
factorial, or Pochhammer symbol, 

{x)j = x{x - 1) ■ ■ ■ {x - j + 1), 

expanding ( polynomial in x we have 



[X 



v=0 



where is the (unsigned) Stirling numbers of the first kind; for instance, 

1 



1' 


= 1, 


j 


J. 




J - 1. 



(j — 1)! and 



0. 



Regarding Stirling numbers of the second kind, denoted {^}, we make use 
of the identity 



k 
j=0 



m 
k 



In particular we note that 

= for m < k, 
Furthermore, 

fe-i 

p=0 



1, and 



m 
k 



k + 1 
k 



(8) 



k + 1 
2 



(9) 



m + 



m 

—y 

+ 1 ^ 

p=0 



171+ 1 
P 



(10) 



where Bp denotes the p*'* Bernoulli number defined by the explicit recurrence 



E 

j=0 



m+ I 
j 



B, = S{m) 



(11) 
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where S{m) is the Kronecker delta function, taking the value 1 for m = 
and zero otherwise. 

For any integer m let ra^mix) denote a polynomial in x of degree no more 
than max(m, 0) whose coefficients are a function of ck, possibly a vector; the 
polynomial ra^mix) is not necessarily the same at each occurrence. Note with 
this notation that ra.,m{x) + „-j(x) = ra,m{x) and if /(a) is a function of 
a not involving x then f{cx.)ra^m{x) = r^^la;). From f lTT]) one easily finds 
that Bq = 1 and Bi = —1/2, and hence ( !T0|) yields 

fe-i 

= ~ 2^" + (12) 

p=0 

Other examples that illustrate this notation are 

fc^ - 2km + + f/m - 2j = r^k,m),2ij) 

and 

fc^ - 2km + m = r^k,m),o{j)- 

Similarly, let ra.^m,n,p{x,y) denote a polynomial in x and y whose coef- 
ficients are a function of a where the highest power of x that occurs is at 
most max(m, 0), the highest power of y that occurs is at most max(n, 0), and 
the highest sum a + 6 in terms of the form x°'y'' is at most max(p, 0); again 
the polynomial ra. jn,n,p{x,y) is not necessarily the same at each occurrence. 
Examples that illustrate this notation are 

k^ + k^ + kY+f + f = n,,MJ,k) 

and 

k^ + mk^f + k^f +f = r„,6,4,6(j,fc) 
An identity that can be obtained directly from Equation 5.114 in [1] is 

n=0 ^ ^ 

for m > 1. We can manipulate this last expression to deduce 

E(-i)"L7i ="("^-^) (14) 

n=0 \ "I" / 
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where u{n) is the unit step function defined to be 1 when n > and is zero 
otherwise. 

Using the well known identity 

and ([T^ we easily find 

/ m — 1 \ , , , 

--6{m-l). (16) 



n=0 ^ 



n 



The following can be found in [3]: 



vn—\ 



, n J n + B T (m + B) 
where, r(m) denotes the gamma function defined by 

t-oo 

Y{z) = / e-^e-'dt, {^z > 0). 
Jo 

We note from [1] that for positive integers m, 

r(m) = (m-l)!, rQ)=v^, r0)=^, (18) 

and 

^ / 1\ (2m- 1)!! ^ 
r m + - = ^ — 



2 y 2" 

where !! denotes the double factorial and when m = we define (—1)!! = 1. 
We can write the last result as 



/ 1\ (2m -1)! ^ , , 



Using ( 1161) and ( 1171) we get 

""^ ^m-l\ n ^, ^, /3r(m)r(/3) 



n=0 ^ 



5(m-l)---^--^M(m-l). (20) 



n J n + l3 r(m + /3) 



Combining f ll7p and f l2U]) we obtain for real numbers a, S, 7, /3 and integer 
m such that 7 7^ and ^ 7^ 0, —1, —2, . . . , — (m — 1), 



5{m-l) + —. ^-[5 - —\uim-l). (21) 



We now derive another result that will be useful. Using (jl2ll and the 
binomial theorem we can write for integers n,p,s > 0, 

k—m 

q{q + ir+^ 

q=l 

k—mpn+s 

= EE 

q=l j=0 
k—m 

= E (^'"^'^' + (^^^ + ^)^'"'"' + r(n,p,.),pn+.-l(g)) 
9=1 

(22) 



^pn+s+l-j 



pn + s + 2 
{{pn + s)(3 - 2m) + 1 - 2m) 
2(pn + s + 1) 

~l~''"(n,m,p,s),p7i+s (^) • 

If we let k = j in this last identity and set m = 1 we get 

^^^^ ' pn + s + 2 2{pn + s + l) ^ ^ 

~^^{n,p,s),pn+sij) ■ 

For the following definition we adopt the empty sum convention that 
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Xq = when b < a. 



q=a 
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Definition 2.1. For nonnegative integers k,j and m, let 

Pk,o{j) = 1, 

k—m j—1 

Pk,m{j) = ^ qPk,m-i{q + 1) - ^qPk,m-i{q + 1) for m>l, 

q=l q=l 

and 

k—m 

P°o(j) = l and PlmU) = J2l^L-i(l+^) form>l. 



q=3 



For example, since Ylq=i 1 ~ ^(^ " for £^11 k >0, for m = 1 we obtai 



am 



The following lemma shows that the identity Pk,i{k) = is a specific instance 
of a more general fact. 

Lemma 2.1. For all nonnegative integers k,j andm, 

Pk,Uj) = forO<k-m + l<j <k, (24) 

and 

Pk,m{j) = Plmij) for 0<j<k. (25) 
Furthermore, for all integers k,m and j satisfying < m, j < k, 

Pk,Uj) = i-^r i^- ' l^^, j + r(.,m),2m-2(j), (26) 

and for integers j and v satisfying < v < j , 

PjA^) = + ^^,2v-2{j)- (27) 



Proof: First note that (IMl) implies (1251) . as the equality holds for < j < 
k — m hy construction. To argue by induction, we have that holds when 
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m = 0, as in this case the premise 0<A; — m + l<j<A;is vacuous. Now 
assume (1241) holds for m > 0. From Definition 12. for k — m < j < k, 

k—m—l j — 1 

Pk,m+iij) = ^ qPk,m{q + 1) - qPk,m{q + i) 

q=l g=l 

= - qPkAq + ^) 

q=k—m 
j 

= - Yl {q-^)Pk,m{q)- 

q=k—m+l 

For k — m + 1 < q < j < k we have Pk.m{q) = by fl24p . Hence, 
Pk,m+i{j) = for k - m < j < k, 



which is fl2^ with m + 1 replacing m, completing the proof of the first two 
claims. 

To prove the rest of the lemma we will show that 



n=0 

1-1 

3 



n ■2n ;,2m— 2n 



i(2n - 3)(_l)-+l^-2n+ipm-2-2n 

2™(m-l-ra)!ra! 

n=0 ^ ^ 

l{2n -2m- l)(_l)nj2npm-l-2n 

2™(m-l-n)!n! 

~l~^m,2'm— 2,2m— 2,2m— 2 

from which both (1261) and ( 1271) will easily follow. 
Now let 

= 2^(771 (29) 

We will use induction on m in our proof. Clearly (128!) is true for m = 0. Now 
assume (l28l) is true for some m — 1 > 0. Then, to evaluate 

k—m j—1 

Pk,m{j) = lPk,m-i{q + 1) - (lPk,m-i{q + 1) (30) 

q=l q=l 
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let 

Pk,mM = Yl 2"^{m-n)\n\ 

n=0 ^ > 

PKm-M - 2^ 2^{m-l-n)\n\ 

n=0 ^ ' 

i(2n - 2m - l)(-l)"j2«A;2'»-i-2« 
P^mM = 2^ 2-(m-l-n)!n! 

n=0 ^ ' 

^m,2m— 2,2m— 2,2m— 2 

so that (!30l) becomes 

k—ra 4 J — 1 4 

^fc,m(j) = X] ? X] Pk,m~l;tiq + 1) " 5 i'fc,m-l;t(g + !)• 

9=1 t=l <j=l t=l 



We will evaluate the terms in fl3T|) separately. We find for t = 1 

k—m j — 1 

,m— 1;1 + 

"""^ r 1^nl,2m-2-2n '^"'^ 
„ -'-^n,m— 1 1 



\2n 



^' l^nL2m-2-2n ■'"^ 
ri=0 g=i 
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Using ([22D and ([22]) with p = 2 and s = 0, ([22D becomes 

k—m j—1 

qPk,m-l;l{q + 1) - ^ qPk,m-l;l{(l + 1) 
q=l q=l 

P"" (-1)'^ A;^'"-! (-l)"(2n(3 - 2m) + 1 - 2m) 

~ 2"^ ^ (m - 1 - ra)!(n + 1)! 2"" ^ (2n + l)(m - 1 - n)!n! 

^1 (_l)n^-2„+2pm-2-2„ (2n - 1 ) j2"+lp™-2-2n 

2'»(m - 1 - n)!(n + 1)! ^ 2"(m - 1 - n)!n!(2n + 1) 



+^m,2m-2,2m-2,2m-2(j) ^) 
m-1 / X 



2™m! 

n=0 



™ 1 /m-A 2(3-2m)ra + l-2m 



n N / 



2™(m-l)!^^ ^ \ n J 2n + 1 

^ ' n=0 ^ ' 

- (-l)"j2nA;2™-2n ^ ^1 ^ _^)n_^-2n+ipm-2~2n _ 1) 

"^2^ 2™(m-n)!n! ^ 2'^(m - 1 - n)!n! (2?2 + l) 

n=l ^ ^ n=0 ^ / V I ; 

~t~^m,2m— 2,2m— 2,2m— 2 

Using dUl) and (El]) with a = 2(3 - 2m), 5 = 1 - 2m, 7 = 2, and /3 = 1 this 
becomes 

k—m j—1 

^ gPfc,m-i;i(g + 1) - X] qPk,m-i;i{q + 1) 

^"(^ - 1) - t;^;:^ TTT w , n ^("^ - 1) - (3 - 2m)5(m - 1) 



2'^m! ^ ' 2'^(m- 1)! I r (m+ I) 

™ (-l)«j2nA;2™-2n ^1 ^ _i )„^-2n+ipm-2-2n - 1) 

+ 2^ 2"'(m-n)\n\ ^ 2™(m-l-n)!n! (2n + l) 

n=l ^ ^ n=0 ^ 7 V ' / 

~t~''^m,2m— 2,2m— 2,2m— 2 
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Using (^Si> and ([11]) we get 

k—m J — 1 

^ qPk,m-l;l{q + 1) - ^ qPk,m-l;l{q + 1) 
9=1 9=1 
/ 1 \n ■2nL2m-2n 

= ^yzlLUl uim-i) 

?i=0 ^ ' 

) u{m - 1) + — -, rrfS - 2m)5(m - 

(2m- 1)! ^ ^ 2'»(m- 1)!^ ^ ^ 

(_l)n_^-2n+1^2m-2-2n - 1) 

- 2-(m-l-n)!n! (2n + l) + "™-2™-2,2„.-2,2™-2(j, fc). 
We next find from fl3Tl) witli t = 2 



k—m j — 1 

^ qPk,m-l;2{q + 1) - ^ qPk,m-l;2{q + 1) 
q=l g=l 

m-2k~m ;^2^ _ 3)(_^^)„+l^^^^ ^)2„+ipm~4-2n 



EE 



3' 



2"^-i(m - 2 - n)\n\ 

n=0 q=l ^ ' 

!^^^ l(2n - 3)(-l)"+^g(g + l)2n+lfc2n^-4-2n 

n=0 (j=l ^ ' 
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Using fl22p and fl23p with p = 2 and s = 1 we get 

k—m j—1 

qPk.m-l-M + 1) - 'iPk,m-l;2{q + 1) 

i^2n-3)(-l)"+ip'"-i 



2^-^{m - 2 - n)\n\{2n + 3) 

|(2n - 3)(_l)n+l^-2n+3^2m-4- 

~ ^ 2'»-i(m-2-n)!n!(2n + 3) 



n=0 

+'"m,2m-2,2m-4,2m-2(j) k) 
1 u2m-l m-2 



2'"-i(m - 2)! M n J 2n + 3 

' |(2n-3)(-l)"+lj2-+3A;2™-4-2n 

2™-i(m-2-n)!n!(2n + 3) 



~^~'^'m,2m— 2,2m— 2,2m— 2 

Using f l2T]) with m = m — 1, a = 2, 6 = —3, 7 = 2, /3 = 3 and then using 
f lTSjl and ([19]) we obtain 

k—m J— 1 

,m— 1:2 + 

q=l q=l 

- ^ ' ^5(^-2) (34) 



(2m -1)! ' ' 2'^-i(m-2)! 

|(2n - 3)(_l)-+l/-+3fc2m-4-2n 

~ ^ 2"^-i(m-2-n)!n!(2n + 3) 

ni,2m-2,2m-2,2m-2(j, 
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Next we find from fl3ip witli t = 3 

k—m j—1 

^ qPk,m^i-M + 1) - X] qPk,rn-l-fi{q + 1) 

q=l q=l 



EE 



3^ 



2™-i(m - 2 - n)!n! 

^ 1^2^ _ 2m + l)(-l)"g(g + l)2n^2m-3-2n 

2-^2^ 2"^-i(m - 2 - n)!n! 

5=1 n=0 ^ 

Using fl22p and (125]) with p = 2 and s = tliis becomes 

k—m j — 1 

qPk,rn-l;z{q + 1) " X] qPk,m-l;z{q + 1) 
9=1 g=l 

" 9--i(^ _ 2 - r2)!n!(2n + 2) + ^^'^^-^l^) 



n=0 



i(2n - 2m + l)(-l)"J2"+2fc2'»-3-2- 
- A. 2— l(m-2-n)!n!(2n + 2) + ^m,2m-3,2™-3,2rn-2(j, fc) 

i^m-i ^2 /^_2\2n + l-2m 



E(-)"C":^ 



2™(m-2)!^^ ' \ n n+l 

|(2n - 2m - l)(-l)"J2"A;2™-l-2- 
+ 2^ 2"^rm-l-nVn' ^ '^m,2m-3,2m-2,2m-2(j, A:). 
n=l ^ '' ' 



Using ( I2TI) witli m replaced bym — l,a = 2, 5 = 1 — 2m, 7 = 1 and /3 = 1 
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this becomes 

k—m 



i-1 



^ qPk,m-i;3{q + 1) - qPk,m-vM + 1) 

9=1 

T(m- i)r(i) 



9=1 



1 l,2m-l 



2™(m- 2)! 



r(m- 1 + 1) 

|(2n - 2m - l)(_l)nj2npm-l-2n 

2™(m-l-n)!n! 

n=0 ^ ^ 

i(-2m-l)A;2"^-i 



(1 - 2m - 2)M(m - 2) + 2(5(r 



2'"(m- 1)! 



U(m - 1) + r„,2m-3,2m-2,2m-2(j, ^) 



- 2m - I)(_l)-j2nfc2m-l-2n 



E 

n=0 



2'"('m- 1 -n)\n\ 



-5(m-2) + 



2™(m-2)! ' ' 2™(m- 1)! 

+'"m,2m-3,2m-2,2m-2 

Next we find from (j3Tll witli t = 4 

k—m j—1 



S{m - 1) 



qPk,m-iAq + 1) - qPk,m-iAq + 1) 

fc— m 

= XI ^ (?^m,2m-4,2m-4,2m-4 (g + 1 , /c) ) 
9=1 

j-1 

— X^ Q' ('"m,2m-4,2m-4,2m-4(Q' + 1, k)) 
9=1 

' m,2m— 2 ,2m-2,2m-4,2m-2 
= ''"m,2m-2,2m-2,2m-2 
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Combining (ESD, dMD, (ES]) and (El]) we get 

k—m j — 1 

^ qPk,m-iiq + 1) - qPk,m-i{q + 1) 

9=1 9=1 
( -\\n ■2nl,2m-2n 

= ^lz£li_^ u(m-l) 

?i=0 ^ ' 

-u(m-l) + — — (3-2m)(5(m- 



(2m-l)! ^ ^ 2'»(m-l)! 

(_l)n_^-2n+1^2m-2-2n - 1) 

- 2-(m-l-n)!n! (2n + l) + "™-2™-2,2„.-2,2™-2(j, fc) 



n=0 



om-l/' _ 1 \lL2m-l lL2m-l 

2 m Ijfc _ _ 

(2m -1)! ^ ' 2^-1(771-2)! ^ ' 

g |(2r7 - 3)(_l)n+lj-2n+3pm-4-2n 



n=0 



2'"-i(m-2-n)!n!(2n + 3) 



~l~^m,2m— 2,2m— 2,2m~2 

"'"^ K2n - 2m - l)(_l)n^-2npm-l-2n 



ra=0 



2™(m - 1 - n)\n\ 



IP^-i ^ ^ i(2m + l)P™-i , 
3 5(m-2) + ^^ p ^— — 5(m-l) 



2'"(m-2)! ^ ^ 2'"(m- 1)! 

+'"m,2m-3,2m-2,2m-2 (j, ^) + ,2m-2,2m-2,2m-2 

Note that for all integers m 

- ^ ^ M(m-1)H -(3-2m)5(m 



(2m -1)! ^ ' 2™(m-l)! 

+ -T-, uim - 2 5- -5 m - 2 

(2m- 1)! ^ ' 2'»-i(m-2)! ^ ' 

^ ^ i(2m+ l)A;2™-i ^ 



2"^(m-2)! ^ ^ 2"^(m-l)! 
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and 



^ 2™(m-l-n)!n! (2ra + l)~^ 2'»-i(m - 2 - n)!n!(2n + 3) 

^1 (_l)n+l 1(2^ _ 3)^-2n+ipm-2-2n 

~ ^ 2™(m-l-n)!n! 

n=0 ^ ^ 



SO noting that m — 1 > in this induction proof ( l37|) becomes 

fc— »n j — 1 

qPk,m-i{q + 1) - X] qPk;m-i{q + 1) 

^ -^'■jnj2n ^2m~2n 



2™(m - n)\n\ 

(2r2-3)(-l)"+ij2"+ifc2 



n=0 

l/o„ nM T \n+l A2n+1 u2m-2-2n 



2™-i(m - 1 - n)\n\ 



n=0 

|(2n - 2m - l)(_l)nj2npm-l-2n 

2'"(m - 1 - n)!n! 

n=0 ^ ' 

~t~^m,2m— 2,2m— 2,2m— 2 



which proves ( |28|) . 

Now we may write as 



Pk,mij) 



l^_Y^mj2m ^_-^^nj2nj^2m-2n 



2™m! ^ 2'™(m-n)!n! 

n=0 ^ ' 

|(2m - 5)(-l)™j2— 1 |(2n - 3)(-l)"+ij2"+iA; 



n=0 

2m-l "^-2 1^^ ^Nj^ T \n+l A2n+1 u2m~2-2n 

2™(m-l)! ' ^ 2™(m-l-n)!n! 



n=0 

|(2r2 - 2m - l)(-l)"j2"A;2"-i-2« 

+ 2'»fm- 1 - n)!n! ^ '^m,2m-2,2m-2,2m-2(j, 



n=0 



-1) ( ) + r(fc,m),2m-2(j) 
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which proves (12^ . Now interchanging k and j and letting m = v 
becomes 

^ (_l)npn^-2.-2n 1 (3^ _ 3) (_ l)n+ipn+l^-2.-2-2n 

~ 2^ 2^(t;-n)!n! ^ 2^(t; - 1 - n)!n! 

n=0 ^ ' n=0 ^ ' 

^ i(2n - 2t; - l)(-l)"pn^-2-i-2" 

+ / , IZr, rn ^ ^ t.,2?;-2,2i,-2,2«-2 (fc, J j 

^-^ Z"[v — 1 — n)\n\ 

n=0 ^ ' 

SO 

J^- - H2^ + 
= + ^"-,2.-2 0) 



which proves fl271) . The proof of the lemma is now complete. □ 
Lemma 2.2. With g'wen m Definition \2.1\ 

{fNN),=J2i-^yPUWNNy-'' forj>0, (38) 

and 

{N~k),^k = Y.(-irPl{k)N=-^-'^ forj-k>0. (39) 

s;=0 

Proof: First note that the validity of fl5Ul) establishes (1551) since with k = 0, 
P°^(0) = P°^(l). Expanding (x — /c)j in x G M for A; G N and j > we obtain 



(x - A;), = ^(-l)X+i-M(A:)x^■-^ 

v=0 

where 

I) 

u<li<l2< - <lv<j m=l 

with any empty product set to 1. Comparison with (1381) shows that it suffices 
to prove 

a,+,_i,,(fc) = P^^ik). (40) 
20 



Equality holds for v = 0. Assuming (HQ]) holds for some f — 1 > 0, we have 

V 

k<h<l2< - <lv<k+j — l m=l 



V 

I'm 



k<q<k+j~v l2,--;lv-q<l2<---<lv<k+j—l m=2 

= E ' E ri'" 

k<q<k+j-v h,...,lv-i:q+l<li<l2<---<lv-l<k+j-l ni=l 
k+j—v 

= X] qak+j^i,v-iiq + I)- 

q=k 

As -Pj°i,(fc) is characterized by the same recursion, the inductive step is com- 
plete, proving fl39|) . □ 
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Proof of Theorem 1.1 

Expanding out the product in ([S]), by (jS]) and we have 

Corr(fc) = e( J](J^-/^)) 

\AeH J 



GCH \AgG 



(N-\G\\ 



GcH \nj 

k /A^-jA 
V n—j J 



3=0 GcH,\G\=j 



E \jJ \n-jJ j n ^k-j 
— jm — 

i=0 \n) 



IN 

k 

\k-3 



where 
and 



fN)'J^^(^fj{-iymjJN), for /^G (0,1) (41) 



In) = f^\fNN)j{N - for J = 0, 1, 2, . . . , fc. (42) 

We also write 

Corr(A;, /^) = 4^ for /^G (0,1). (43) 
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Using ( 155]) and (15^ . and recalling that = when f < then 

becomes 

/ j \ fk-j 

\v=0 / \i=0 

k / k—r \ 

= fNj2[ Y.(-^yP'A^^fN''{-^)'-''-'Plk-.-rij))N^ 
r=0 \v=j—r / 

= E E fN''pim-^r-"pik-v-r{j)N^ 

r=0 v=j—r 

k k—v 

= E fNP'vi^) E i-^f^'pL^^-rm"^ 

v=j—k r=j—v 

= T.fNPli^) E (-l)'-^nV.-.o•)iv^ 

v=0 r=j—v 

since -P°^(l) = for f < 0. 

For a function represented as the power series = J^JL-oo^jf''^ ^'^^ 
k & Z let fi{f; k) denote the coefficient of f'' in that is, 

Making the change of variable m = A; — r in the first sum below, and 
again using that PkmU) = for m < 0, for < < A; we obtain 



k—v 
r=j—v 

V 

= pU^) E i-^rpL-.m'-"- 

•m=k—j+v 
k—j+v 

= p"A^) E i-^rPk,m-.ij)N'- 

m=v 

k—j+v 

= E Kk,J,fN;-v:m) (44) 
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where 

X{k,j,f^; -v:m) = (-l)"i^?,(l)P°„0')iV'"™- 
Now, by dUD and OS]), 

k 



a{kjN] k-v) 



j=0 

k k—j+v 



j=0 m=v 

k+v k—m+v 



m,=v j=0 
k+v k—m+v 



iyX{k,j,h;-v) 

{-iyX{k,jjN; -V : m) 

m=u jr=0 ^ ^ 
in=v j=0 ^ ^ 



Tk~m 



Tk—m 



where in the last two equahties we invoke Lemma \2A] to apply Pkm~vU) = 
for J > k-m + v, and = Pfc°„(j) = Pfc,„-.(j) for < j < fc. 

At this point in our proof we will consider k even and k odd separately. 
For k even it will be convenient to write and (B7I) . respectively, as 



and 



= (-1) 



m 3 



2m 



2™m! 



+ T{fc,m),2m-l(j), 



■2v 



+ rv,2v-lij)- 
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Substituting these last two expressions into (HSj) we get 

k+v k 



a{kjN] k-v) 



m=v j=0 



^ ^'^m—vj2(m—v) 

k+v k 

;-i)'-EE (^K-i) 



j+m 



X 



./ \. ~^ ^(k,m,v),2m—l 

'.[m — v)l 



AT' 



k—m 



From 0431) and (HSl) we obtain 



j^rk/2 (^i'^jN;k-v) 



jSjk/2 



{N)k 

k+v k 

D^EE 

m=v j=0 



X 



-^^m—vj2m 

2'^v\(m - v)\ 



,m,v),2m—l 



N' 



k—m 



which becomes 



N''^^Coir{k, /at; k-v) 

j^k/2 (^_l^k~v k+v 



1 

2™(m - v)\ ^ \j 



j A2m 



+ Y1 [■]i-'^yrik,m,v),2m-lij) 
3=0 ^-^^ 



N 



k—m 
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Using ([H]), (H7I) becomes 

N''/^CoTT{kjN;k-v) 



^''■'(-^r-Cf 1 (48) 



k+v k , , \ 

+ E E ( J (-l)^'^(^,".,^),2n.-l(j) N'^-'^- 
m=v j=0 ^ / 

Since (A^)fc is of degree k, then when m > k/2, 

J^k/2 ATk pTk/2 

hm ——N^-"" = hm — = 

and when m < k/2, from (jH]) and we have that 

^ (''^ {-lyf = {-l)''k\ {f } = for all p < 2m - 1, 

and therefore, by linearity, 

P j (-l)^?^(fc,m,«),2m-i(j) = for m < /c/2. 



i=o 

and since for /c even 

'2m 
k 



for m < k/2, 



then upon letting — )■ oo, f l48p is zero for f > k/2 and for f < k/2 
[Hi) becomes 

lim A^'=/2qq^^(^^ j^. ^ _ ^) (49) 



N^oc{N)k v\ 1 2"\m-v)r ' \k (> 



k—m 



When m = k/2. 



T\Tk/2 yyfc /Yfc/2 /yfc /yA:/2 

lim A^''"'" = lim = lim = 1 
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Therefore, using 



becomes 



lim N''/^ ■CoTT{kjN;k-v) 



i-iy 



k\ fk/2 

V 



2^/2 (I)! 

{-lyEz 



klk/2 



Now, by our previous notation, Corr(A;) = Corr(A;, /^r), so 
lim N''/^ ■ CoTT{k) = lim A^^/^ • Corr(A:, f^) 



N^oo 



lim N''/^ V Corr(fc, k - v)f 



k/2 



v=0 



v=0 



= [f{f-l)]'/'EZ' 

which proves the theorem for k even. 

For k odd we have from (ESI), (EZD and (USD that 



a{k, fN]k - v) 

k+v k 



:-i)^EE (-1) 



2^w! 



m=v j=0 



-1) 



+'"(fc,m,D),2(m-t))-2(j) 



2™-''(m -f)! 



AT' 



k—m 
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which becomes 



k+v k 



m=v j=0 



-ly- 



X 



-ly 



-V /^i^™ + I [(^ - v){2m -2v-5) - v{2v + 



2m- 1 



+?^(fc,m,-!;),2m-2(j) 



AT' 



k—m 



From ( 143|) and ( 146|) we obtain 



fc+?; A: 



9?(-i)^eeGi(-i)- 



X 



-ly 



~v ff"" + i [(^ - v){2m -2v-5)- v{2v + 



2m- 1 



+'^(fe,m,t;),2m-2(j) 



N' 



k—m 



^ ^ i=o ^-^ 



E (-1) 



X ( j + 3 [(^ - ^)(2m - 2t; - 5) - v{2v + 1)]^' 



2 m- 1 



j=0 



N 



k—m 



28 



Using ([H]), (13UI) becomes 
iV('^+i)/2Corr(A;, fN]k-v 



2'"(m -f)! 



-1)'=A;! 



2m 

A; 



(51) 



^ ^ |[(m-t;)(2m-2t;-5)-t;(2t;+l)] ^ 

+ 2^^ ^wz^ — — (-1) ^' 



2'"(m - i;) 



fc+D k 



2m -1 
k 



fk\ 



k—m 



m=v j=0 

Since (A^)fe is of degree k, then when m > (fc + l)/2, 

A^{A;+l)/2 Tyfc /v(fc+l)/2 

hm — -— AT'^— = hm — = 

and when m < {k + l)/2, from and we have that 

k 



-ly/ = (-l)^A;! {^} = for all p < 2m - 2, 
and therefore, by linearity, 

-l)^'^(fe,m,^),2m-2(j) = for m < (A; + l)/2, 
and since for k odd 

I 9m 1 

= for m < (A; + l)/2. 



j=0 



and 



2m 

A; 



2m- 1 

k 



for m < + l)/2. 
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then upon letting N ^ oo flSip is zero for v > {k + l)/2 and for v < {k + 
becomes 



lim N^''+'^^/'^CoTT{k, Jn; k-v) 



= lim 

(fc+l)/2 1 
m=(fc+l)/2 

When m={k + l)/2, 

]\l{k+l)/2 



(fc+l)/2 

E 

m={k+l)/2 



1 



2'"(m - 1;)! 



-l)''k\ 



2m 
k 



[(m - v){2m -2v-5)- v{2v + 1)] 
2'"(m - 1;)! 



-1)'=A;! 



2m- 1 

k 



lim 



k—m 



lim 



jyfc jY(fc+l)/2 



Therefore. fl52l) becomes 



jyfc ^(fc+l)/2 



lim Ar(^+i)/2 . Corr(fc, - 



-1) 



k—v 



1 



2(fc+i)/2 -t;)! 



-1)'=A;! 



k + 1 
k 



+ 



I [(^ -v) {k + l-2v -5) - v{2v + 1 



2('=+i)/2(Mi_^)! 

Using (|9]) this last result becomes 

hm Ar('=+i)/2 . Corr(fc, fN]k-v 
\{-lYk\ 



{-lfk\ 



y!2(fc+i)/2 (Mi _^)! _ 
'k + 1 



k + 1 
2 



+ 



-v]{k-2v-A) -v{2v + l) 



\{-iy{k-l){k + l-v)k\ 
t;!2(fc+i)/2 -t;)! 
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We may write this last result as 



lim iV(*^+^)/2 ■ Corr(A;, fjv, k - v) 



^ l{-iY{k-i){k + i-v)k\{f) 

2(fe+l)/2 (fe±l)! 

l{-iy{k-i){k + i-v)Ez''+^ 

~ k+1 \v 

Hence, 

lim ArC^+i)/^ . Corr(A;) = lim A^e^+i)/^ . Corr(A;, /jv) 

N-i-oo N-^oo 



= lim TVC^+i)/^ V Corr(fc, h; k - v)!^' 



t)=0 



k + 1 \ V 



k + 



2 



A; + 1 df 

= [/(/ - 1)]^'=-^)/'(2/ - - 

which proves the theorem for k odd. 

The proof of the theorem is now complete. □ 
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